PAIR PRODUCTION 
IN A TIME DEPENDENT MAGNETIC FIELD 



Giorgio CALUCCI* 

Dipartimento di Fisica Teorica dell'Universitd di Trieste, Trieste, I 3401 4 Italy 

INFN, Sezione di Trieste, Italy 

Abstract 

The production of electron-positron pairs in a time-dependent mag- 
netic field is estimated in the hypotheses that the magnetic field is uniform 
over large distances with respect to the pair localization and it is so strong 
that the spacing of the Landau levels is larger than the rest mass of the 
particles. This calculation is presented since it has been suggested that 
extremely intense and varying magnetic fields may be found around some 
astrophysical objects. 
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1. Statement of the problem 

Standard Q.E.D. clearly states that a static electric field of sufficient 
strength and extension can create pairs of opposite charges through a tun- 
nel effect [1,2], but the same kind of production is not foreseen for magnetic 
fields in the same conditions, unless one would take into account the pro- 
duction of monopole pairs. The situation changes if the external fields vary 
in time and the effect of the field variation on the pair production rates 
has been discussed in detail in connection with the possibility of detecting 
these effects in experiments in optics [3], where the time dependence of 
the field is explicitly taken into account, with an essential negative an- 
swer, at least for the real pair production.* It may be possible that pairs 
are created in the presence of very large field strengths with time varia- 
tion which is not too small; since it has been proposed that situations of 
this kind could take place in some phases of neutron star evolution[5,6], 
a preliminary investigation on the pair production in a very strong and 
non constant magnetic field is presented here. No detailed quantitative 
calculation is performed, but some order of magnitude estimates are given 
in order to show that the model calculation that will be presented can be 
well suited for the conditions referred to above. 

These conditions correspond to a very intense magnetic, over 10 10 T, 
which remains constant over distances which are large in comparison with 
the dimensions involved in the electron production, during the time of 
variation the spatial uniformity is conserved, in this first investigation 
only a change of strength at constant direction is considered. 

The rapidity of change cannot be too small in order for the production 
process not to be entirely negligible, but since it involves macroscopic 
objects it will be slow enough to allow the adiabatic treatment of the 
production process [7]. For the ratio between the magnetic and the rest 
energy, defined as 9ft = (2heB) / (m e c) 2 , we obtain for the chosen values of 
B, 9ft > 4.5, which allows some computational simplifications that will be 
used. 

2. Determination of the production rate 

The production rate is calculated here in the following external field 
configuration: the magnetic field grows at fixed direction, this could hap- 
pen locally, in the case of a collapsing neutron star, when a compression 
of the magnetic field takes place [6]. 



* The effect of virtual pairs in external fields is a different issue, which is under 
investigation also experimentally [4], but it will not be discussed here. 
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In order to fix notations and conventions some well known features of 
the relativistic electron in static uniform magnetic field are rewritten here. 
The Dirac equation* Hifj = [a • (p + eA) + flm]^ = et/j is written as: 

( m a-(p + eA)\ ( <f>\ = ( <f>\ 

\a-(p + el) -m J \x J \X J 

and the system is solved for </> when e > and for \ when e < . 

When the field becomes time dependent, the choice of the vector po- 
tential for the magnetic field becomes physically relevant, different choices 
yielding the same magnetic field give rise to different electric fields. So 
a choice is unavoidable: when the magnetic field is taken as constant in 
direction and varying in strength a cylindrical configuration, defined by 
A = — if A S, will be studied. Decomposing the variables into transverse 
and longitudinal components (the longitudinal direction is called z), the 
equation for is 

\vl + Pt + \e 2 r 2 T B 2 + eB(L z + a z )](j) = [e 2 - m 2 ]</> 

The variables p z and a z are diagonalized with eigenvalues k and s 
with the position 

= ^= e ifc2 ^, s , 

we obtain 

[Pt + h 2 r 2 T B 2 + eB(L z + s)]0 fejS = (e 2 - k 2 - m 2 )(j) k}S . 

For future use the solution of the remaining equation is expressed in term 
of the harmonic-oscillator operators. Following ref. [8] one can define: 

a>v = (Pv/b — ibv)/V2 with v = x, y b 2 = \eB 

ad = (a x - ia y )/V2 , a g = (a x + ia y )/V2 , N d = a\a d , N g = a g a g 

and express in this way the operators pt , ty , L z . The final result is : 

e 2 =m 2 +k 2 + eB(2N d + 1 + s) . 

The operators N d , N g are standard number operators, we respectively 
call n , p their eigenvalues and remark two known results relevant for the 

here e > so the charge of the electron is — e. 



next developments: the quantum number p has no influence on the value 
of the energy; the lowest Landau level, defined by n = and s = — 1, has 
an energy independent of the strength of B. 

The second order equation for x is the same as the equation for 0; 
the normalized solutions of the original equation are: 

j>+=jV( t 6 ) > i>-=jV(~^A (6) 

V = a ■ (p + eA) , iu = |e| , N = y / (w + m)/2w (6') 

The second quantized electronic field is decomposed in terms of these 
solutions 

vI/ >>.><••., , * t = E c M,7+^<7 ( 7 ) 

so the standard representation of the (second quantized) Hamiltonian is 
produced: 

H = f d 3 r^(r) [a • (p + eA) + (3m] *(r) = wj [c] Cj - djd]} . (8) 

The index j embodies all the relevant quantum numbers j = (k,n,p,s), 
of course when they are read as positron quantum numbers, the sign of k 
and s must be reversed. 

Now the variation of the Hamiltonian with time is given writing 
B — Bk and B = B a + Bt and defining h = dH/dB and a± = ^(a x ±ia y ). 
In terms of the harmonic-oscillator operators, eq.(4), it results: 

h = — \ea A f • k = i— [— a + a,d + ct-a d — a-a g + a + a g ] , (9) 

The derivative of the second quantized Hamiltonian is then: 

|^ = / <Pr&h9 = 
dB J 

3,3' 

(1<») 

When we are interested in the transition from the vacuum to the electron- 
positron pair configuration only the second of the four addenda is in- 
teresting and will therefore be calculated explicitly; all the addenda are 
important for higher order corrections. 
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Before going on with the calculation of the matrix element, the fol- 
lowing observation (ref. [8]), is useful: the mean value of r\ is found to 
be 

< n, p\r T \n, p >= (n + p + l)/b 2 

In semiclassical states not having sharp eigenvalues the mean value be- 
comes * 

« n, p\r\\n, p {yfp - \fn) 2 jb 2 . 

In both cases one observes that n cannot grow too much because it costs 
energy eq.(5), while p can grow without cost; so if the space at disposal of 
the system is large in most configurations we shall find 

< n, p\r\\n 1 p >~ p/b 2 . 

It is, therefore, convenient to define p = b 2 R 2 and R may be interpreted 
as the mean position; the sum over p which at the end occurs, is converted 
into an integration according to J2 P ~ ^ 2 / dR 2 . The region of integration 
will certainly be large, as compared with the microscopic lengths, so the 
terms containing R will be the most important and they will be calculated 
using the approximate representation a g ~ bRe lS , a g ~ bRe~ ld . The 
conclusion is that the dominant contribution for macroscopic extension of 
the magnetic field may be extracted after reducing the expression of h to 
the much simpler form 

h Q = ±ieR(e- [9 a + - e ie a_) . (11) 

The matrix element corresponding to the transition from the vacuum 
to a two-particle state is represented as J d 3 rip^_ jh ip-j> = \ieRNjNj> M, 
two quantum numbers are equal in j and j', they are k and p. Using the 
relation: 

V = a z k + bV2[a + a d + cr_a^] 
the expression for M is found to be 

M =\(a + e- w - a-e [9 ) - ^ -x 

l. [Wj + m)[Wj> + m) 

[-k 2 (a+e- w - a-e w ) + V2kb(a z a\e~ ie - a z a d e i0 )+ ( 12 ) 
2b 2 (-a+a d a d e w + a -a\a\e- ie )] } . 

* this procedure is similar to the treatment of a condensed Bose system, where 
one puts a ~ a) ~ vaJa. 
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From this form it appears that the two produced particles cannot have 
wholly equal quantum numbers, in particular it is impossible that both are 
produced in the ground state. There are eight relevant matrix elements 
which are listed as follows. 





=< 


-,n\M\+, n >= 


— € 


T 2 


=< 


+,n\M\-,n >= 


e~ 


T 3 


=< 


+,n|M|+,n+ 1 


>= 


T 4 


=< 


+,n + l\M\+,n 


>= 


T 5 


=< 


-,n\M\-,n+ 1 


>= 


T 6 


=< 


+,n+l\M\+,n 


>= 


T 7 


=< 


+,n\M\-,n+2 


>= 


T 8 


=< 


-,n + 2\M\+,n 


>= 



iO 1.2 AO \ 



= -V2kbe i9 Vn + ![(«;_,„ + m)(w- >n+1 + to)] 1 
y/2kbe~ ie Vn + T[(w- jn + m)(w- n+l + m )\ 
2b 2 e i6 y/(n + l)(n + 2)[(w-, n+2 + m)(w+, n + m)]" 1 



-2b 2 e-' ld y/(n + 1) (n + 2) [{w + , n + m)(w- ,n+2 + Tn)\ 

(13) 

Further simplification is gained for strong magnetic fields, in fact in 
this case the excited Landau levels are very far from the fundamental one 
so the production of an electron-positron pair in their ground state would 
strongly be favoured, but it has just been seen that this transition is not 
allowed and at most one particle may be found in the ground state of its 
transverse motion. We may neglect the mass term in comparison with the 
magnetic energy and write: 

w +j7l = a/46 2 (?i + 1) + k 2 , w-, n = a/46 2 ?! + k 2 if n > 
W-fi = Wk = V m 2 + k 2 . 

The dominant terms for large b of the matrix elements are independent of 
b and they are explicitly: 

Ti = - T 2 * w -e ie 

^ _ (13') 

T 5 = - T* w -e ld k/(w k + to) if n = T 5 = -T* -> if n > 

f 7 = - T 8 * w ±e ie . 

The same approximation is introduced into the normalization factors and 
for them it results J\f l/\/2, but for jV-,0 = a/ (wfc + m)/2wk- The 
amplitude for the pair production is written in terms of the standard adi- 
abatic expression [7], with the following specifications: only the transition 
from the vacuum to the two particle state is considered; the reduction of 
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amplitude of the vacuum state is not taken into account; the initial instant 
of the process is not pushed to — oo, but it is some finite time, t = 0, in 
which the field B(0) = B Q is already large. This is expressed writing the 
adiabatic projection coefficient [7] 7 J;J / as: 



7j 



m R C a <j,j'\dH/dB\ > . f 

,f(t) = -Bj dr — 7 exp[-iy AEdt \ : 



with AE = Wj + Wj', eq (14). 

Prom now on we refer to the matrix elements Ti, T2, in the other cases 
slight modifications are required. 

It is convenient to separately treat the n = and the n > contribu- 
tions. In this last case, in which the sum over all n will be performed, we 
shall set n n + 1 and in so doing 

£ AEdt' = [(2enB(t) + k 2 f/ 2 - (2enB a + k 2 f /2 ] ; 

the second term gives rise to an irrelevant constant phase, in the calculation 
of \"ij,j'(t) \ = \eRB\J\ the relevant integral is 

r* dr ^ r — ^L (2enS(r) + A .2 )3 /2 

Jo 2^/2enB(r) + k 2 L 3enB 
A convenient form for this expression is 

J =— \f(S ) ^2enB + k 2 + f(S(t)y/2enB(t) + k 2 )] 
2enB 

/oo 
exp[-iqy 3 ] dy 

2 , (15) 

S = ^(2enB + k 2 ) 3/2 

3enB 

S(t)=-^(2enB(t) + k 2 f 2 . 
SenB 

The parameters S are very large in all the interesting situations. This is 
better appreciated if we insert the powers of c, % that are needed to come 
back to standard units, for simplicity , when k = 

s= 4cvm /2 

3B 

then from the general asymptotic expansion of /, 



the first term is sufficient * and we obtain: 

J = -ii[(2enB + fc 2 ) -1 exp[iS ] - (2enB(t) + k 2 )' 1 exp[iS(t)}] . 

Moreover, since S(t) is very large, the term containing cos[<S — which 
appears in |J7"| 2 is rapidly varying around and it will be dropped in the 
following. This us allows to write: 

|l71 ~ 16(2enS + k 2 ) 2 + 16(2enB(t) + k 2 ) 2 ^ 

The sum over the longitudinal quantum number is performed with the 
usual transition to the continuum variable — > J dk, the sum over 
the transverse quantum number remains discrete and the final result is 



The symbol £ represents the Riemann's z-function [9]. 

For the estimate of the n = term the difference between the two 
addenda (+) and (— ) is important and we obtain: 

J AEdt' = kr + ^i- [(2enB(t) + k 2 ) 3/2 - (2enB Q + k 2 f /2 ] = $ ; 

The subsequent integral is more complicated: 



J ° = J ¥ e 

but has the same qualitative features as in the previous case, $ is large 
and we consider the asymptotic expansion: 

$2 $4 

It may be verified that $/$ 2 is of the order 1/S. Going on as in the 
previous case, we keep only the first term in the expansion, we drop the 
rapidly oscillating term and we obtain: 

\Jo\ 2 ~ [k + \/2eB + k 2 } " 4 + [k + y/2eB(t) + k 2 } " 4 (18) 



* it is possible to give a close expression for /, in terms of the irregular confluent 
hypergeometric function U [9], although this is not very useful. 



The sum over the longitudinal quantum number yields: 



2Z 

15tT 



+ 



(2eS ) 3 /2 (2eB(t)f/ 2 



Now we must simply collect the terms, i.e. the factor ^{eRB) 2 , relating 
\J | 2 to I7I 2 and the integration \eB J dR 2 , which performs the sum over 
the remaining quantum number p, the final result is: 



1 



1 1 



B(t) + 1 



16 V ^L256 SV " / ^ ' \^\ ZR ' M(?l2B Y B ll 2 ' y /B^j. 

The production rate in a cylinder of height Z and radius Rm, i n a uniform 
and very strong magnetic field, varying linearly with time is given approx- 
imately by dl/dt. If we evaluate dl/dt we obtain a cubic dependence on 
B, built up by an original linear dependence of the transition amplitude, 
which is standard in the adiabatic processes and has to be squared, and by 
the dependence on B which is mainly related with the fact that the pro- 
duced state extends over a transverse region becoming smaller and smaller 
as the magnetic field grows. 

While the dependence on the longitudinal variable is trivial the de- 
pendence on the transverse one, i.e. the fourth power, is less obvious and 
it may be understood if we remember that the configuration is not trans- 
lationally invariant in the transverse dimension. In fact, although the 
magnetic field is uniform, the electric field £ grows linearly with R and 
this puts an evident limitation to the configurations one can consider: in 
fact too large values of R such that S = \BR > B cannot fit in this treat- 
ment. It is pleasant to find a numerical confirmation of the intuitive idea 
that the n = configuration must be important, in fact \J Q \ > The 
rest of the factors strongly depend on the details of the chosen configura- 
tion. In particular there must be a further and not very large factor which 
takes into account the other transition-matrix elements Tj , i > 2 . 



3. Conclusions and outlook 



A calculation of the pair production in a time-dependent magnetic 
field has been presented, in a simple situation where the already strong 
field is still growing in a fixed direction, B > 0. This process is suggested by 
the consideration of an imploding neutron star where the magnetic field is 
locally compressed [6]. The choice of a very strong field has a less dramatic 
effect on the production rate than one could have expected at first sight. 
In fact, to the growth of the state density in x-space there corresponds a 
growth in the distance among the energy levels and therefore a decrease 
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of density in momentum space. A stronger production rate could take 
place if both produced particles could be found in the ground state, but 
the dynamics does not allow this process. We may note, however, that as 
the distance of the Landau level grows there is an increasing number of 
longitudinal levels between two transverse levels, the convergence in the 
longitudinal momentum is always rapid, as dk/k 4 , so the dependence on k 
which has been neglected in the terms Tj cannot play any important role. 
A more formal and general observation can be done: the total angular 
momentum J = J d 3 r f A [£ A B] is zero in the considered volume and field 
configuration while an e~ e + pair in the lowest energy level has no orbital 
momentum and opposite spins, so the total angular momentum is 1, the 
transition is therefore impossible. 

Another interesting configuration could be given by a magnetic field 
which changes its direction, like e.g. that one could find around a rotating 
neutron star [5]. The adiabatic formalism is certainly still useful, a first 
look at the problem shows that in this case some difference in the treatment 
is needed since the adiabatic states are, of course, time dependent, but the 
energy levels remain constant. Anyhow a more detailed investigation is 
required in order to obtain a definite answer. 

Much more complicated is the possibility that an intense magnetic 
field may act together with the QCD vacuum leading to an enhanced 
production of quarks and so, finally, of hadrons. Similar effects have in 
fact been proposed [10], but any quantitative estimate is very uncertain. 
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